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Abstract 

We discuss the far-field spatio-temporal cross-correlations of 
waves multiple-scattered in a turbid medium in which is em- 
bedded a hidden heterogeneous region (inclusion) character- 
ized by a distinct scatterer dynamics (as compared to the 
rest of the medium). We show that the spatio-temporal cor- 
relation is affected by the inclusion which suggests a new 
method of imaging in turbid media. Our results allow quali- 
tative interpretation in terms of diffraction theory: the cross- 
correlation of scattered waves behaves similarly to the inten- 
sity of a wave diffracted by an aperture. 



A considerable progress has been made during the re- 
cent years in tlie understanding of wave transport in dis- 
ordered media.u Very similar phenomena are shown to ex- 
ist in multiple scattering of electrons and classical waves 
(e.g., light) under particular circumstances J3 Some of the 
concepts developed first theoretically, and then studied 
in model experiments, are now very close to practical ap- 
plications. One of the important fields where the physics 
of multiple-scattered waves is currently finding its ap- 
plications is the (medical) imaging of disordered, turbid 
media.u The light waves scattered inside a turbid medium 
(e.g., human tissue) carry information on the properties 
of the medium. The information can be considered as 
being "encoded" in the statistics of the waves. Analy- 
sis of the latter statistics allows one to reconstruct (or 
"image") the scattering medium. 

A simplified version of a typical geometry considered 
in connection with imaging problems is shown in Fig. |^. 
A slab of turbid medium occupies the space between the 
planes z = and z = L, and some region (a cylinder- 
shaped inclusion) inside the slab is assumed to have 
somewhat different properties as compared to the sur- 
rounding medium. If "different properties" means differ- 
ent scattering /i^ and/or absorption ^.a coefficients, one 
can image the inclusion by measuring the spatial distri- 
bution of the average intensity /(|Cl|= {E{r,t)E*{r,t)) 
of transmitted (or reflected) wave.QB Here E{r, t) is the 
amplitude of scattered wave at spatial position r at time 



t. If /i^ and /ia are constant throughout the medium, 
and the contrast between the inclusion and surround- 
ing medium is provided by the scatterer dynamics (dif- 
ferent types and/or intensities of scatterer motion inside 
and outside thg inclusion) , the methods of diffusing- wavg 
spectroscopyQ'Ll can be applied to visualize the inclusion.^ 
In the latter case, one measures the time autocorrela- 
tion function Ci (r, r) — {E{r, t)E* (r, t + r)) of scattered 
wave field at multiple positions r, which allows visualiza- 
tion of the inclusion. 
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^Published in Waves and Imaging Through Complex Media, 
edited by P. Sebbah (Kluwer Academic Publishers, Dordrecht, 
2001). 
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FIG. 1. A cylindrical inclusion of height 2A = Z2 — zi ^ £ 
and radius a 3> ^ (^ is a photon transport mean free path) is 
embedded at zo = (zi + Z2)/2 inside a turbid slab of width 
L and surface area A = W^ IF > L, IF > a. ka and k^. 
denote the wave vectors of incident waves, while kj, and kj, — 
the wave vectors of transmitted waves. The scatterers in the 
medium undergo Brownian motion with diffusion coefficients 
Din (inside the inclusion) and Dout (outside the inclusion). 

In the present contribution, we propose to use the 
spatio-temporal cross-correlation function Ci(r, Ar,T) — 
{E{r,t)E* {r -\- Ar,t -\- t)) for the purpose of imaging 
in turbid media. Since the time autocorrelation func- 
tion Ci(r,T) carries more information about the turbid 
medium than the average intensity /(r), we suggest that 
the information contents of the spatio-temporal correla- 
tion function Ci(r, Ar, r) should be even more rich. If 
the points r and r -I- Ar are taken far enough from the 
medium (in the far-field of scattered wave), Ci(r, Ar, r) 
is equivalent to the angular-temporal correlation func- 
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tion {E{kb,t)E*m,t + T)) [where E{k,t) is the spatial 
Fourier transform of E{r',t) with r' taken at the plane 
where the scattered waves leave the medium] . 

We start with a macroscopically homogeneous turbid 
medium (no inclusion) , and assume that scatterers in the 
medium undergo Brownian motion with a diffusion co- 
efficient D. In addition, wo assume a weak-scattering 
limit k£ ^ I (where £ — l/zij is the photon transport 
mean free path), and neglect the absorption of light in 
the medium {^a = 0) . As depicted in Fig. [l], a plane wave 
is incident upon a turbid slab at time t with a wave vec- 
tor ka . The transmitted wave leaves the slab with a wave 
vector kb. Similarly, for an incident wave with a wave 
vector kjj at time t + r, the transmitted wave has a wave 
vector k'f,. Assuming unit amplitudes of incident waves, 
we calculate the correlation function of transmitted fields 
CiiK,kb;K,Kr) = {EiK,kt,t)E*{K^lt + T)} us- 
ing the standard diagrammatic techniques Jj following the 
general calculation scheme developed in Ref. In the 
leading order in a small parameter l/k£, we obtain: 



d^Ri d^Ra 



X exp (-jAqaRi + iAqbR2) 
X P({Ri,n,{R2,L-n,r), 



(1) 



where q's denote projections of k's onto the plane z = 
const, Aqa = qa - qjj, Aqt qf, - q|,, r = {R, z} with 
R being a two-dimensional vector perpendicular to the 
0-axis, and we assume the first and the last scattering 
events to occur at z = ^ and z ~ L — £, respectively. If 
L ^ £ and |ri — ^ £, the reduced ladder propagator 
P entering into Eq. (m), obeys the diffusion equation: 



a 



(r)] P(ri,r2,T) 



^3 



S{ri - ra), 



(2) 



where a^ir) = 3t/{2to£^) with tq = {APD)-^. The 
solution of Eq. (||) with Dirichlet boundary conditions at 
z = and z PI L (F = if zi = 0, L or 22 = 0, L) is 
readily foundO 

Po(ri,r2,T) 



smh[f3a{L - z>)] sinh(/3aZ<) 
Pa smh{(3aL) 
X exp [i{Ri - R2)p] . (3) 

Here = + oi^(T), z> — max{zi,Z2}, z< — 
min{zi,Z2}, and the subscript "0" of Pq denotes macro- 
scopically homogeneous case. Inserting Eq. (^) into Eq. 
(0), we get 



Aqb pj= the time autocorrelation function of transmitted 
lightQ is recovered. The Kroneqker delta symbol in Eq. 
(^ describes the memory effect .113 

Now we turn to the case of macroscopically heteroge- 
neous medium, assuming that the scatterer diffusion co- 
efficient Din inside a cylindrical region depicted in Fig. |l| 
is not the same as Dout in the surrounding medium (while 
£ is assumed to be constant throughout the whole sam- 
ple). The correlation function of transmitted waves can 
be again described by Eqs. (Q), (||) but with a^ir) — 
ag(r) -I- aKr) = 3T/{2Tin£'^) inside the inclusion and 
a'^ir) = a^(r) = 3t / {2Touti'^) outside it. Here nn.out = 
i4PD,n,outr' and aliT) = 3t / i2£^)[l / - l/ront]- 
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FIG. 2. Normalized angular-temporal correlation of a wave 
transmitted through a macroscopically homogeneous (no in- 
clusion) slab of width L = 101 for three different time delays 
r/ro — 0,0.01,0.1. This correlation function vanishes iden- 
ticalljj-^or Aqa 7^ Aqt, which corresponds to the memory 
effect .y 

Assuming |ai(T)| <C UqIt), we can write an approxi- 
mate solution of Eq. (^ as a sum of Pq corresponding to 
the macroscopically homogeneous medium with tq = Tout 



[see Eq. J 
inclusion:! 



and Pi which describes the influence of 



Pi({Ri,£},{R2,L-n,r)c. 

/ d^r Po({Ri,n,r,r) 
X Poir,{R2, L-£},t) 



144^3 A 
X exp(ipRi — isR2), 



d2pd2sFT(p,s,r) 



(5) 



(4) 



where the integration of the second line is taken over the 
volume of inclusion, and Ft is a form factor: 



'1 (Aqa,Aq.,r)^JA,„.Aq.^ ^^^^.^^^(°^^ 
with = Aql + a^T). For a2(T) ^,0, Eq. (|) reduces ^ , ^ ^ 2Ji(a|p-s|) 



to the angular correlation function,0 while for Aq^ ~ 



a IP 
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(MiM sinh(/3,L) sinh(/3bL) 
r sinh [iPa - l3b)A] 



sinh \(3a 



cosh [/3a (i - Zo) + PbZo] 

cosh [Pa{L~ Zo) - (3bZo] 



(6) 



Here Ji is the Bessel function of the first order, 



«§(t), P'b 



+ a5(T), A = (z2 - 2;i)/2, and zq 



(zi + z2)/2 (see Fig. |lD. 

Inserting P = Pq + Pi into Eq. we obtain the 
angular-temporal correlation function corresponding to 
the macroscopically heterogeneous medium as a sum of 



two contributions: Ci 
by Eq. (|) with ^(t) = 

C['\Aqa,Aqb,T) 



CI 



(0) 



C'[^\ where c[^^ is given 



q;o(t), and 



X —jFTiA(la,Aqb,T) 



(7) 
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FIG. 3. Normalized angular-temporal cross-correlation 
function C[^^ corresponding to the setup of Fig. |l| We assume 



Aqa II Aqt, for this plot. In contrast to Cj"'' (see Fig. C} 
is not necessarily zero for Aqa 7^ Aq;,. 
r / 0. 



^(1) 



C[^^ ^ only for 



Equation (^ is the main result of the paper. We 
now compare the C^"'' correlation function, correspond- 
ing to a macroscopically homogeneous slab [Eq. (^], 
and the c[^^ correlation function [Eq. (0)], originating 
from the presence of a dynamically heterogeneous re- 
gion (inclusion). As follows from Eq. (^, C^"'' van- 
ishes identically if Aq^ T^^qt, which is a manifesta- 
tion of the memory effectO If Aq^ — Aqf,, c[^'' de- 
cays to zero for Aqa > 1/L (see Fig. ||). In contrast, 
c[^^ correlation is not necessarily zero for Aqa ^ Aqf, 
(see Fig. The memory effect is still present for the 

C^^'' correlation function, as it is peaked near Aqa — 
Aqb due to Ji(a [Aqa - Aqb|)/(a [Aqa - Aq^l) term in 
FT(Aqa, Aqb, r). The memory effect for C^^'' is consid- 
erably less sharp than for C^"' , as one can see from Fig. ||. 



Let us consider the simplest and practically important 
case of a single incident plane wave (Aqa = 0). For a 
macroscopically homogeneous slab, the angular-temporal 
cross-correlation vanishes if Aq^ ^ 0, i.e. the waves scat- 
tered in different directions are uncorrelated.^ If a het- 
erogeneous region is embedded inside the slab, the 
term appears and correlation between the waves scat- 
tered along different directions is not necessarily zero. 
The C[^^ ter m is plotted in Fig. ^ for three different radii 
a of the inclusion (solid lines). As is seen from the figure, 
the correlation range can be estimated as Aqb ~ 1/a. It 
is worthwhile to note that the correlation between the 
waves scattered along different directions (Aq^ 7^ 0), in- 
troduced by the inclusion, exists only for t ^ 0. If r = 0, 



CI 



(1) 



and the correlation function is given by Cl 



(0) 



which is 


identically zero 


for Aqf, ^ Aqa = 


0. 
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FIG. 4. Normalized angular-temporal correlation functions 
of transmitted waves for a single plane wave (Aqa — 0) in- 
cident upon a slab with a cylindrical inclusion inside (solid 
lines). Dashed lines show the (normalized) angular distri- 
bution of the wave field £ diffracted by a circular aperture 
of radius b {b = 5A and b — 20A with A being the wave- 
length). For a = 20£ ^ A and 6 — 20 A, correlation and 
diffraction curves are remarkably close, suggesting that cor- 
relation is "diffracted" by the inclusion. 

Some qualitative insight into the behavior of the c[^^ 
correlation can be gained by comparing Eq. (^ with the 
angular distribution of the wave field £ (wave number 
K — 2tt/K) diffracted by ^a, circular aperture of radius b 
(see Fig. ^, dashed hnes):l 



£{K±_) cx b 



,MK^b) 



Kj_b 



(8) 



^In reality, correlation persists as long as A^i, < l/W, and 
our result m) corresponds to the limit W ~* 00. 
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where K± is the projection of K onto the plane of the 
aperture. In the case of a ^ A, which corresponds to a 
"pill-shaped" inclusion, the correlation functiongiven by 
our Eq. and the diffraction pattern of Eq. (|| are re- 
markably close (see, e.g., the curves corresponding to a = 
20£ and b = 20A in Fig. In this case, one can explain 
the appearance of correlation between the waves scat- 
tered in different directions in a macroscopically heterp* 
geneous medium using the classical diffraction theory,t^ 
and assuming A = £. As a consequence, some theorems 
known for diffraction of waves (e.g., the Babinet's prin- 
ciple) , apply directly to the angular-temporal correlation 
function of light transmitted through a macroscopically 
heterogeneous turbid medium. It is worthwhile to note 
that this holds for any shape of inclusion, provided that 
the transverse extent of inclusion is significantly greater 
than its extent along the z-axis (a 3> A in our notation). 
For a ^ A, the quantitative agreement between Eqs. iM) 
and (^) is absent, although their overall behavior is sim- 
ilar (see, e.g., the curves corresponding to a = 5^ and 
= 5A in Fig. |). 

Up to now, our analysis has been devoted to the 
correlation functions of transmitted light. In experi- 
ments, however, it could be more convenient to work 
with diffusely reflected waves. Calculation of the angular- 
temporal correlation function of reflected waves is per- 
formed similarly to that of transmitted ones. If |qa + | , 
iQa + QbL |qa + Qbl, iQa + I > ^c Can ignore the 
time-reversal symmetry and obtain: 



^(Aqa, Aqfc,T) =(5Aq„.Aq,-^ 

^ smh[pa{L-e)]smh{pai) 



X FR(Aqa, Aqf,,r), 



(9) 



(10) 



FR(p,s,r) 



2Ji (a |p — s|) 



a|p-s| 
sinh(/3a^) sinh(/3b£) 



iPaiW) smh{(3aL) sinh(/3fcL) 
f sinh [{(3a + (3b) A] 

sinh [{pa - (3b) A] 



{(3a ~ (3b)A 



■cosh[(/3a ~ (3b){L ~ zq)] 



(11) 



The cj"'' correlation function given by Eq, 
the result of Ref. 



The ^ term [Eq. 

.(1) 
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|9|) reduces to 
for r = 0, (3aL oo and (3ai 0. 
(p^j has the same qualitative fea- 
tures as the C\^' correlation of transmitted light [Eq. 
(^]. In reflection, however, low-order scattering events 
become important for Aq^, Ag^ ^ l/£, and thus our re- 
sults (p|)~(pT|) make sense only for Aqa(,Aqbi ^ 1. If 
|qa + q&l, jqi + q|,l, |qa + c^\, or + qf,| are of order 



of or smaller than the time-reversal symmetry of the 
problem cannot be ignored any more. This significantly 
complicates calculations even-in the case of macroscopi- 
cally homogeneous medium .113 

In conclusion, we have calculated and discussed the 
angular-temporal cross-correlation functions of waves 
scattered in a turbid, dynamically heterogeneous 
medium. Our analysis demonstrates that the considered 
correlation functions can be used to image a hidden dy- 
namic inclusion embedded in an otherwise homogeneous 
medium. Comparison of our results with diffraction pat- 
terns obtained for a wave diffracted by an aperture sug- 
gests that the angular-temporal correlation function can 
be considered as being "diffracted" by inclusion. Such an 
interpretation is particularly successful for "pill-shaped" 
inclusions which are much more extended in the trans- 
verse directions (i.e., in the directions parallel to the sur- 
faces of the slab where they are embedded) than in the 
longitudinal one. 
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